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THEOjRETICAL SUPERSONIC NAVE iseAG OF UHTAPERED SWEFTBACK 
AND EECTMCaiAE WINGS AT ZERO LIFT 
By Sidney M. Harmon • 


SUMMARY 


A theoretical investigation of the supersonic wave drag at zero 
lift of a series of untapered sweptback wings having thin symmetrical 
biconvex parabolic-arc sections has been iresented in NAQA TH No. 1319* 
The investigation has been extended to include Mach numbers which 
bring the Mach line behind the wing leading edge and also to include 
wings of rectangular plan form. The results are presented in a uiiified 
form so that a single chart permits the direct determination of the 
wave drag for this family of wings over an extensive range of sweep- 
back angle, Mach number, aspect ratio, and thickness ratio. The 
results obtained for the total wave drag of the sweptbaok wings are 
applicable to the same family of wings having a coi’responding degree 
of sweepf orward . 

When the Mach line lies behind the wing leading edge, the wave- 
drag coefficients of the sweptback and rectangular wings are shown to 
reach maximum values at certain limiting aspect ratios and remain 
constant for all aspect ratios greater than these limiting values . 

The limiting aspect ratio is equal to — =— — for the 

' cot A - 1 - 1 

sweptback wing and to ■■ for the rectangular wing, where A 

\/m2 - 1 

is the angle of sweepback and M is the Mach number. The variation 
of wing wave-drag coefficient with Mach number over the complete 
range of supersonic Mach number is shown to beccaae less pronounced 
with decreasing aspect ratio. It is also shown that sweepback obtained 
by rotating the wing panels rearward can give appreciable reductions 
in wing wave-drag coefficient at all supersonic speeds. 

INTRODUCTION 

Recent developments in airfoil theory for supersonic speeds 
(references 1 to 3 ) indicate pronounced effects of sweepback and 
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aspect ratio on the drag. In reference 3 , a method haaed on thin— 
airfoil theory for a frictlonleas fluid (reference 2 ) waa applied to 
calculate the supersonic wave drag at zero lift for a series of wings 
having thin symmetrical biconvex parabolic~arc sections with untapered 
plan forms and various angles of sweepback and aspect ratios. The 
results in reference 3 ^ however, are limited to those cases in which 
the Mach line lies ahead of the wing leading edge, or to a range of 
Mach number from 1.0 to the value equal to the secant of the angle 
of sweepback, The term "Mach line" as used herein refers to the 
Mach wave that originates at the loading edge of the center section, 
unless specified otherwise. 

The present paper extends the calculations of reference 3 for 
the same series of wings in order to present wave— drag results for 
cases in which the Mach line lies behind the wing leading edge. These 
data are obtained for wings of ractanguldr and swoptback plan forms 
and cover an extensive range of Mach number beginning with the value 
at which the Mach line’ .coincides with the wing leading edge . Jn order 
to present a more ‘c camp late picture, the results of reference 3 that 
cover the lower range of Mach number ar.^' reproduced herein together 
with additional cui’ves ccaaputed from foniulas given in reference 3 . 

The results of the entire investigation are presented in a unified 
form similar to thatngiven in. reference 3 so that t'he wave drag for 
this .family of wiiigs may be determined directly from a single chart 
over an extensive range of sweepback angle, Mach number, aspect ratio, 
and thickness ratio. In the Tenth Annual Wright Bros. Memorial 
Lecture given on December 17, 1946, Dr. von Ka'rraan indicated that at 
zero lift the total wave drag for a sweptforward wing is identical 
with that obtained for the same wing having a corresponding amount of 
sweepback. The results of the present investigation for the total 
drag’ of sweptback wings, therefore, are applicable to the same family 
of wings having a corresponding amount of. sweepf orward . The, distribu- 
tions of. section drag, however, will differ in the two cases . Although 
the calculations have been made for. the biconvex parabolic--arc profile, 
the data may be applied to Ihdicate corresponding results for profiles 
similar to the biconvex parabolic— arc profile. 


SYMBOLS . 


X, y, z coordinates of mutually perpendicular system of axes 
c chord of airfoil section, measured in flight direction 

t/c thickness ratio of section, measured in flight direction 

A angle of- sweep, degrees 
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m = cot -A' ;Vr’ ■ Ij,. ■ — 

- •-'- ^ -y _ '■ 

h wing semiapan. measured along y-axls, semichords except in 

appendix A-. .' ■ -.r-y'-r^. %f -* / . ■■ ■ . •■ 


K 

.A-. , 
S ■ 

I ft 


peranister indicating span\Siaa position equal- to y/m,' sekichords 
■■aspect ratio ‘ ■ 


wing -area ■ ' ‘ 

. . , • . . > 1 ’ . . . ’ , •■ , ^ • • - 

slenderness ratio, ratio of wing semispan measured along 
midchord' station to maximum thickness of center sedtion 


V velocity in flight direction ’ , 

.u X— component of disturhanoe velocity, positive in flight 

, direction 

u distiirhance velocity caused hy source, line with reversal in 

' .- ■■■ sign of m ■ .... , 

I source factor required to maintain a given wedge angle 

^ slope of airfoil surface 

E.P. real part of complex expression 

M. Mach; number . • ■ ' ' ■■•; - ' ■ ' 

"0 = /m2' 1 ‘ . ’ ' 

coordinate measured along y— axis which is shifted to tip 
section, semichoi’ds 

y-jj . coordinate maasured along y-axis which is. shifted .to 
opposite tip section, semichords 

Cdoo^.- , .section -wave-drag coefficient without tip effect 

c^ section wave-drag coefficient including tip effect 

Ac^ increment in section ^wave-drag coefficient caused by wing “tips . 

Acdj -increment in ’ section wave-drag coefficient caused by'Ving tip 

located on same half .of wing as . section 



4 


MCA UN No, 1449 


Ao^ increment in section wave,~drag coefficient on one Ving panel 
caused ty tip of opposite wing panel 

wing wave-drag coefficient without tip effect 

,Cp , wing wave-drag coefficient Including tip effect 

ACp increment in wave-drag coefficient caused by tips, complete wing 

ACj)^_ incrementr-in wave-drag coefficient on one wing panel caused 

by tip of opposite wing panel, complete wing 

I , Ti coordinates which replace x and y, respectively, used to 
indicate origin cf source line 

v' V in transformed coordinate system of reference 2 

primed values of A,, y, y^, y-jj, h, z, and m indicate transformation 
involving multiplication by factor j3. 

Subscript notations for u and indicate the origin of so^Irce line 
in terms of coordinates x and y, respectively, 

ANALYSIS 


; 


Basic data. - The apaiysls is based on thin-alrfoil theory for 
small pressure disturbances relative to the ambient preeeures. The 
axes used are the mutually perpendicular x, y, z system in which 
the x-axls is taken in the direction of flight positive to the. roar, 
the y-axls is along the span positive to-the right, and the z-axls 
is positive .upwards. The symbols used to designate the wing— plan— 
form parameters are shown in figure 1. The analysis is made- for 
untapored sweptback and I’eotangular wings of biconvex parabolic— arc 
profile at zero lift. The wing' is considered to be cut off in a 
direction parallel to the- directlc«n of flight. The Mach numbers 
considered in this analysis correspond to those for which the Mach 
lino lies behind -the wing . leading edge, that is, for Mach nimibers 
greater than the secant of tho angle of sweepback (m> 


Theory .— The present analysis corresponds essentially to that 
given' in reference 3 where m^ By' following the analysis of 


reference 3 , the section wave-drag coefficient for the ; symmetrical 
biconvex parabolic-arc profile at the spanwlse station y is 
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■where u refers to either surface of the airfoil. 

The evalmtion of eq.uation (l) involves the determination of 
the disturbance velocity u at each point on the wing. An appropri- 
ate system 'of semi— inf inite source and sink lines which represent 
the wings considered herein is given in appendix A. Appendix A also 
gives the u-e3cpressions foy these sowce lines .and their regions of 
influence on the wing as determined by the Mach cones from their 
extremities. These regions of influence for the individual source 
lines are illustrated in figure 2 and are given in table I as the 
limits of the variables of integration for x along the chord and 
for y along the span. 

Figure- 2 is given in order to. illustrate typical Mach lines for 
Mach numbers in which the Mach wave frcsm the nose, lies behind the wing 

leading edge . ^m > Inasmuch as the wing cu-fc— off is represented 

by reversed semi— infinite source— line distributions (appendix A), the 
tip Mach cones in figure 2 for the various aspect ratios show the 

extent of the region of the 'tip. effect. For A> — SS — , the tip’ 

' mp + 1 

effect on' the wing is influenced only by the adjacent tip . If 

A < both wing tips influence each wing panel, 
p 

Formulas for section wave-drag coefficients. — The formulas for_ 

the section wave-drag ooefficlents for sweptback and rectanguiar 
wingSj which result .frcm the integration of the u-expressions in 
equation (l), are presented in- appendix B. These formulas give 
expressions for section wave-drag coefficient without the tip effect 
c^^ and also the expressions for the increments in section waye-drag 

coefficient caused by the wing tips Ac^. 

Wave-drag coefficients for complete wing .— The formiilas fcr the 
wave-drag coefficients for untapered sweptbac’k and rectangular wings 
of biconvex parabolic-arc profile are given in appendix B. 


BESULTS AND DISCUSSION 


Variatio n of section wave-drag coefficient along span.— Figures 3 
and 4 show the variation of section wave-drag coefficien-fc c^ along 
the span for the wings ojf and 0 ° eweepback, respectively. 
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The results are shown in figure 3 (a) for the sweptback wing 
of infinite aspect ratio, that is without the tip effect, at several 
Mach numbers. The lowest Mach number l.Ul represents a case in which 
the Mach line is ooincident-uvith the wing loading edge. In this case, 
as noted in reference 3 , the wing has a very high da*ag and the section 
wave-drag coefficient increases in' the outboard direction and 
approaches infinity at an infinite, distance from the wing center. 

As noted later in this section the assumptions of the linearized theory 
may be violated, at a Mach number bf.'l.kl except for very thin wings, 

-t ■ ’ • ' ' ■ 

For the higher Mach number of 1.50^ '^he ordinates are reduced 
in magnitude and the spanwise variation of cj^ . is markedly flatter . 

The section wave-drag coefficient without tip effect c^ has a- 
minimum value at the center section, incroasen in. the . outboard direction, 
and- then approaches a constant find to . value , This constant value of 
0 ^ results frcm the fact that if the effect of the wipg -tips is • 

neglected, the flow is two-dimensional over the parte of the wing 
ahead of the Mach wave from the vertex.- The cmponent parallel to 
the wing leading edge- has no effect here; that I 0 , this region is 
influenced only by the component of the velocity 'normal ,tp the wing 
loading edge (reference 4). The prese^n’es, therefore, • are exactly 
those that would be oompu'ted by the Apkeret theory of- linearized two- 
dimensional: supereonlc flow. by use of the normal velocity ccsaporient 
(reference 5)* Outboard of the point where the Mach vaye frem the 
vortex intersects the wing trailing edge, the flow^ip entirely two- 
dimensional. In this region, therefore, the section wave-drag coef- 
ficient based on all parameters raeasurod normal -to the wing leading 
edge iias the constant Ackeret value for an. infinite re’ctai^ular wing. 

The value is: . , 



(See appendix B , equation (-B4 ) . ) 


where the subscript n indicates that, p, and t/c are 

lasasured normal to the wing leading edge. The' section wave^rag 
measfired normal to the wing leading edge is obtained frpm this section 
wave-drag coefficient, and the cemponent of this fj?rce in the flight 
direction gives the true eectlpn drag. As a result there is obtained 




cj . cos3 
H»n 



3^m2p9 „ 1 
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When this value for is ccsnpared with the Ackeret result for 

CO 

an infinite rectangular wing having the saioe thickness ratio in the 
flight direction, in the region of twc^imensiohal ’fIow“is 

found to increase because of aweepback by the factor .. 

\/m2p2 „ 1 ■ 

The highest Mach number 2.24 represents a case in which the 
wing is well outside the Mach co ^ iq . At this Mach number, c,i 

'**CQ 

approaches the two-dlmensionai value at a comparatively smll 
distance outboard of the center section. 


In figure 3(b), the spanwise distributions of c^^ are shown 
for the 45° sweptback wing at a Mach number of 1.8 for various 
aspect ratios'. The as]^ct ratios in this figure v^ere selected' so 
that they represent each of the different types of t-lp effects. 
(See fig. 2.) 


A comparison of the c^^^-diabributlons illustrated in figure 3(a) 

with those given In figure 2{«) of reference 3 shows an important 
difference between the conditions where the Meoh line lies ahead of 

the wing leading edge, and where the Mach line lies behind 

the leading edg, (n > i). K n < i. the section w,ve-drag coef- ' 

flclent without tip effect decreases to zero at some point 

along the span, then becomes negative, and finally approaches zero 

asymptotically at infinity; whereas if m > does not decrease 

to zero with increasing values of y but approaches a constant 
positive value. The contributions of the adjacent tip effect to the 
shape pf, the • o^-dlstributlOn are similar' for both 

m>.^ and m < This may be seen by taking the difference of 

'the ■ c^^-curves for A =«> and A = 5 in figure 3(b) and comparing the 
resulting distribution with given in figiire 2(b) of roference 3. 


The increment in section wave-dreg coefficient on one w.'.ng panel 

caused by the tip off the op^slte wlrjg panel Aca , however, tends to 
} . ' ' XI 

be shifted inboard when the Mach number is increased from the range 


m < to m 5* (Note Mach lines from opfosite tip in fig. 2(a) 
of this paper and fig* 2(b) of reference 3-) This com- 

parison indicates, therefore, that in general for a given sweep- 
back angle ahd comparatively high aspect ratios in which Acct.j.j is 

A > j or very sma3.1,. an increase in the Mach number 


zero 
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which hrings the Mach line hahind the leading edge of the wing 
appears to reeult^ln a shift of the center of pressure of the wave 
drag in the outboard direction. 


The results are shown in figure 4(a) for the rectangular wing 
of infinite aspect ratio at several Mach nujnbers. These results 
correspond to the Ackeret theoryj which shows a constant section 
wave-drag coefficient along the span. 


The results are given in figure 4(b) for the rectangular wing 
at a Mach number of 1.25 for several aspect ratios. The aspect 
ratios in this figure, as in figure 3 (b), were selected so as to 
represent each of the different types of tip effects. (See fig. 2.) 


JIffect of tips on wing wave-drag coefflclentr- The present 

1 , - 1 
analysis for ^ indicated, as noted in reference 3 for m < p 

that the integrated value of ACj^^ over the wing is zero if the 


2 m 


Inasmuch as 


aspect ratio is eq.ua! to or greater than 

om mp -i- 1 

ACpj is zero for A « — — (see fig, 2(a).), the total increment 


mp + 1 


in wave-drag contributed by the tip is zero if A 


For the rectangular wing, m »oo and the total increment in 


Cj) caused by the tips is zero if the aspect ratio 
appendix B, equation (Bll).) In the range A ^ i. 


A > (See 
therefore, the 


wave-drag coefficient for the rectangular wing Is independent of 
aspect ratio and is equal to the Ackeret result for a two-dimensional 
wing, If the aspect ratio for the rectangular wing is less than 
l/p, the increment in Op- caused by the tips is foimd to be negative. 


Generalized ciurves for wing wave-drag coefficient .— Figures 5 
to 8 present generalized curves for determining the wing wave— drag 
coefficient over an extensive range of avoopback angle, Mach number, 
aspect ratio, and thickness ratio. The results are given in figures 5 
and 6 for the sweptbaok wings. As noted, previously in the "XNTEOBUCJTIOlf , 
the data in these figures are applicable to the same family of wings 
having a corresponding amount of sweepf oi-vard . The results for the 
rectangular wing are presented in figures 7 and 8 , 

The data in figures 5 '4o 8 apply specifically to uptapered wings 
at zero lift with biconvex parabolic-arc profiles and the wing tips 
cut off in the direction of flight. The results, however, may 
be applied to indicate approximate results for profiles similar to 
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tihe '''bi.ediiT 0 af parato3.i6^rd''p]f'<>^ Beoaufe. ctf ’.the ^se;^ of 

the present linearized thedi'y'j the 'fe suits presehted herein for Mach 
nmhers in the vicinitjif ,qf esQ. v-A, .or .at-.-veryohilgh' Mach ni^hers. 

'are dn. general q^uestiona'ble, except tccc' small 'dr vanishingly email 
thickness , ratios^ inasmuch as ■tjhe.. cad-cuiated. pressure disturhances 
ovhr the Trtng at such Mach ’nimib'ers t.end to he’, large .re^-ative to ,th© 
amhierct pressures; For applicatitjne" to very small aspect . ratios, ■ 
the present theory may r.eq.uir^sbmo’-iflodifioations. ’... . • ■ . i- 


The results in figures 5 6 for the sweptback wing. are given 

in terms of the wave-drag-coefficient parameter the Mach 

mmiber parameter cot A yM“= — 1, and the aspect ratio parameter- 
A tan A, As shown in reference ’ 3 , these parameters group the 
variables Cj), A, A, M, and t/c in a unified form and thereby permit 
the direct determination of .the drag coefficient frcm a sin^e chart 
over an ejctensive range of sweephack angle, Mach number, aspect ratio, 
and. thickness ‘ratio. For a 10— percent— thick sweptback iWir^, the 

drag— coefficient parameter in figures 5 ‘Q^dd 6 becomes simply Cp, 


the aspect-ratio parameter, A, -and the Mach numbed parameter, ‘ /m^ — 1. 
The results in figures 5 and. 6 may be applied to any sveepback angle 
oo'vering a range of asp ect ratio frcm 0 to 10 cot A and a 'Tange of 

Mach number frcm I to /l -i- 49 tan^ a. ■ 

The results in figu res ^(a ) and 5(b) refer to Mach humbers 

corresponding to cot A /m^ — 1 equal to or less ..than 1; that is, 
where the Mach line lies ahead of the wing leading edge. These . 
data represent results which were obtained directly from figure 5 of 
reference 3 with several aOlditiona.l curves which ■were calculated 
frcm the formulae given in reference 3 . The results in figure s 5(c) 

and 5(d) refer to Mach numbers corresponding to cot A '{m- — 1 equal 
to or greater than 1 ,* that is, where the Mach- line lies behind the 
wing leading edge. These data were calculated frcm the forniulas 
given in appendix B of the present paper. The data in figure' 6 were 
obtained by cross— plotting the results in figure 5 ; therefore, the 
range' of Mach number and aspect ratio in figure' 6 is .the same as 
that ^ven in figure 5 . ' 

Effect • of aspect ratio aiid Mach number on wing wave— drag Coef- 
ficient for sVreptback wings .— If the Mach line is we3J. ahead 6 f the 
wirig leading edge, increasing the aspect ratio 3^ the range ' A ^’cpt A 
reduces Cx) (see, fig., 5 (a)) » however,, for all aspect ratios where • 
the Mach line approaches the wi ng leading ..edge arid ’also' for A ^ oot A 

for Mach numbers . cot A — 1 » 1, Cp is reduce^. .wiVh decreasing 
aspect ratio. - w’ \ 
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If the Mach Ijne lies behind the wing leading edge, the wing wave- 
drag coefficient reaches a maxinum value which is constant for all 

aspect ratios gi’eater than a certain limiting value (A ^ 

mp 1/ 

Figure 6 Indicates that In general the variation of Cjj over 
the coiiiplete range pf Mach numbei* becomes lees pronounced as the 
aspect ratio is decreased. For A ^ cot A the variation of C^j 
with M is ccmparativ03.y sauall. 

Effect of aspect ratio and liach number on wing wave-drag coef- 
ficient for rectangular wliig'o .— iFlgures' 7 and. B indicate the effects 
of aspect ratio and Mach number on Cjj for rectarigular wings. 

Cn 

Figure 7 shows the variation of — k with aspect ratio for 

i“(l) 

constant Mach numbers'. The wing wave-drag coefficient reaches a 

maximum value at an aspect ratio equal to l/i/^^ — 1 and remains 
constant for higher values of aspect ratio. 

_ Ct) 

Figure 8 shows the variation of with Mach number M 

100(1)2 

for constant aspect" ratios. When the aspect ratio is equal to or 
greater than 1, Cp is independent of aspect ratio for Mach numbers 
equal to or greater than l.i(l. The curves in figure 8 indicate that 
the variation of Cp over the ccaaplote pai^ge of Mach number becamas 
less pronounced as the aspect ratio is reduced; thus, the same trend 
noted previously for the sweptback wings. is shown. . 

Effect of sweepbeck angle on wing wave-drag coefficient .— In 
order to study the effect of sweepback on the wing wave drag, the 
method of obtaining the sweepback must be ooncldered. In the present" 
investigation, wing wave-drag results have bean pbtained for two 
different methods of increasing, the sweepback angle, and the data are 
presented in figures 9 and 10. ' . 

In the first method of— obtaining sweepback, ' the data for which 
is given in figure 9j the sweepback angle, is increased by rotating 
the wing rearward about a vertical axis at the midpoint of the center 
section. The root and tip sections of phe rotated wing are then 
modified so that- they 'boccme parallel to the fli^t direction in order 
to conform to the family of ' s-weptback wings considered.. here in. For 
this method, the aspect raplq and. the .thiciaiQBs ratio in the flight" 
direction .are reduced with increasing, sweepback,' but the wing area 
and slenderness ratio are maintained constant. The slenderness ratios 
in figure 9' are based on a thickness ratio of 'O.'IO meusitred in a 
direction normal to the wing leading edge. The thicJgaess ratio t/c 
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measured in the flight direction, therefore, varies with sweephack 
as cos A oz* is eq,iial to 0,1 cos A. The aspect ratio is reduced 
with sweephack hy the factor cos^a. The_ aspect ratio A is related 
to the slendetness ratio l/t hy the following formula; . 

" A = 0-2— cos^ A 
t 

The results in figure 9 indicate that if the wing is ■ sweptbaok hy 
rotation on the hasis of. o octant _ Z/t, sweephack can give appreciable 
reductions in the •wave— drag doefficien't at all supersonic speeds. 

For the second inethod of ob’tainlng sweephack, the da-ba for which 
is shown in figure 10, the sweephack angle is increased hy sliding the 
sections rearward. For. this method, the aspect ratio, the thickness 
ratio in the flight direction, and the wing ai’ea are maintained 
constant. The results in figure 10 indicate that if the ^^ing is swept- 
hack hy sliding the sections rearward, sweephack can give appreciable 
reductions in wave— drag coefficient only at I-lach numbers for which 
■the forward Mach line is well ahead of the wing leading edge. At 
Mach nimhers, however, for which 'the forward Mach line approaches the 

wing leading edge or is behind it (cot A \/m~ — 1 ^ 0,95)# wing sweep- 
hack obtained with constant A and t/o inoreases the wave-drag 
coefficient. 


COWCLUDII'fG EEMAEKS 


The theoretical investigation of the supersonic •wave drag of 
untapered swepthack wings at zero lift has been extended to include 
Mach numbers which bring the Mach line behind the wing leading edge 
eind to include wings of rectangular plan form. The wing sections 
investigated were biconvex, composed of two parabolic arcs, and the 
wing tips were considered to be cut off in the direction of flight. 
’The following conclusions have been drawn: 

Mach line ahead of the wing leading edge; 

1. If the Mach line is well ahead of the’ wing leading edge, 
increasing the aspect ratio in the range of aspect ratio greater than 
or approximately eq.-aal to the cotangent of -the' angle of sweephack 
reduces the wing wave-drag coefficient, 

2. For aspect ratios less than approximately the cotangent of 
the angle of sweephack and for all aspect ratios where the Kiach line 
approaches the wing leading edge, the wii;g wave-drag coefficient 
decreases with decreasing aspect ratio. 



.12 


NACA 1’N No.. l41t-9 


■Mach liae •'behind "the wing- leading edgfe: 

3 . The wave’-drag coefficient of -the sweptback wing 'reaches a 
maximum: value at an aspect ratio for which the Mach line from the 
leading edge of the center section intersects the trailing edge at 
the tip; this maximum coefficient remains constant for higher values 
of aspect ratio. 

4. The wave-drag coefficient of the rectangular wing reaches 
a maximum value at an aspect 'ratio for which the Mach line from the 
leading edge of the tip section intersects the trailing edge of the 
tip on the opposite wing panel; 'this maximum coefficient remains 
constant at higher values of aspect ratio. 

5 • Tile wave— drag' coeff.iclent“of all- wings for all aspect 
ratios decreases with increasing Mach number. 

■Complete range of Mach number; 

6 . With a comparatively high aspect ratlOj an increase in the 
Mach number which moves the Mach line behind the wing leading 'edge 
appears ;to move the center of- pressure- of the wave drag in the out- 
board direction. 

7- The variation of the wing wave-drag coefficient with Mach 
number over the complete range of- supersonic Mach number becomes 
less pronounced as the aspect ratio is decreased, 

8 . Sweepback obtained by rotating the wing panels rearward can 

give appreciable reductions in the wave-drag co_efficient at all 
supersonic- speeds . : 

9 . Sweepback obtained by sliding each section rearward, can give 
appreciable reductions in wave-drag coefficient only when the Mach 
line is well ahead of the wing leading edge. When, the Maph line 
approaches the wing leading edge or is behind it, sweepback obtained 
by sliding each section rearward increases the wing wave-drag 
coefficient . 


Langley Memorial Aeronautical Laboratory 

National Advisory Ccmmittee for Aeronautics 
Langley Field, '7a.', July 10, 194? 
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APPENDIX A . . : 

FOEMULAS FOR u-lXPEESSIOWS AND INTEGRATION LIMITS FOR INTEGRAND 
IN EQUATION (1) FOR UOTAPERED SWEPTBACE AND RE^AJteULAE 
WINGS OF BICONVEX PARABOLIC-ARC PROFILE AT Z3iEO 
LIFT = cot .4 = i 
SweptbacE Wings 


The desired integrand- u ih ©^[uation (l) is dete3?mlned in a 
manner similar to that described in. appendix A of reference 3. The 
sweptback wing of desired profile shape and plan form is ‘built up 
by superposition of the solutions obtained for a semi— inf ini te 
oblique wedge. On the basis of the linearized theory, reference 2 
derives a solution representing an oblique semi-infinite (sweptback) 
source line making the angle, of .swoepback A with the y-axis . The 
solution utilized for the' pressure field or for the disturbance 
velocity is ' . 


u, 


0,0 


= R.P. I cosh~^ 


X — mj3^y - 

P(y - mx[ 


(Al) 


where the subscript notation indicates that the source line starts 
at, the origin of coordinates (x = 0,-y = 0). Equation (a1) Is 
shown in reference 2 to satisfy the boundary condition for a thin 

oblique wedge making the half-angle in the transformed ..... 

coordinate s.ystem of refer ence 2 (y’ = yp, = zP)^ where 

^ ^ sotirco factor I in terms of the 

physical coordinate system required to maintain the desired wedge 
is shown in reference 3 to be represented by 


j _ Y m .. dz 

/l - n^p2 ^(0,0) 


(A2) 


If the aouroe line is swept ahead ,of the Mach lino, mp > 1, 
then substitution of equation (A2) in equatibn (Al) gives 
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u. 


0>0 


R.P. 



cosh' 


-1 X - 


y - HZ 


Vm 

as l■ll■l<^l rt i^^ll^ 





cos~^ 

Bjy - 



(A3) 


The swepthach wing of biconvex parabolic-arc profile is 
represented by the followijig systeia of eleaiontary semi-dLnf inite 
oblique source lines i Source lines of eq.ual strength are placed 
along the leading and trailing edges beginning at the center section 
In conjunction with a constant distribution of sinlr lines along 
the chord also beginning at the center section. At the tip, where 
the wing is cut off, reversed semi— infinite soiirce and sink lines 
are distributed so as to cancel exactly the effect of those origi- 
nating at the center section in the entire region of space outboard 
of the tip. 


For m> i, the sources which are considered to originate at 
0 ■ ' • . 
the wing trailing edge havcf no effect on the wing because their Mach 
linea' are behind the wing. (See fig. 2(a) . ) J?ho diaturbance 
velocity u is then expressed in the following form (reference 3> 
appendix A, eq.uatlon (Al) ) ; 


^ “ ’\),0 + ~ H“o,o " D ^ 0,0 


~ '^h/m,h '^h/mj-h ■’' ^'^/m,h ■*' ^h/m,~h 


where the subscript notation indicates the origin of the source 
line. The bars over u refer to the source llhes caused by the 
opposite wing panel; that is, H Indicetes a soiirce lii>^ with a 
reversal in the sign of m. " 

In eq.uation (a 4), the u-expressions are given by the real 
parts of -the following, expressions.:, . . 


u (x, y) = 

n/a2p2 


Vm 






-...Mr."”! “ 5 ■“ 

H|3^(y - n) 

(A5) 

; . , P| X - Tl 

- m(x - i)\ 

.ea-2. ^ + 

mp2(y _ n) 

(A6) 

■ Pjy - n 

.+ h('x - |)j 
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where repreeeuts the origin of the elementary source lines, 

JTor the hiconvex parabolic-arc profilsj 



The symbol i in equation (At) refers to an integration 

which represents the influence of a unifom distributlc^ of source 
lines along the chord of the biconvex profile beginning at the 
position IjT). This integration is expressed as (reference 3 , 
appendix equation (A2) ) : ■ 


7) - 


'x-pjy-ii{ 


^ ejy-rri * 


_ J ooeh' 


dx*^ 


1 


'i„- l X - I 

^[y - Tij 


1 

^ m(x - 1) 


X — t" „„ 

3(y - n) 

*v 

|/m?p2 __ 1 

y - T} 


1 - 



L 


y - t 

* 


'> (A7> 


where I* is the variab3.e of integration representing the t -coordinate 
of the origin of each source line in the distribution of source lines. 
For the biconvex profile. 


d^z 

dx~ 




Equation (A 7 ) is expressed as a function of S-Z - L * that is, 

y - T} • 




y) = ~(y 11 ) f 


-Z-i. 


-(y ~ h) 


Then 
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Rectangular Wings 

■ ■ ’’v'-r ■ -• . 

For the rectangular wirig-, '• bT = 00 and' the following, expressions 
result from substituting m =oo in equation (A5 ), (a 6), and (A7) ; 


i> n np dx 

jt — 

cos“ 

1 P(y - ti)l 

X ^ 1 J 

^ y) = — ~ -cos ^ 

dx ■ X'-'- 

jz.nl 
- & 

■ ... 

iu (x, y) = -Z(y . 

D 

- n) 

d^z j 

X — 1*' 

dx2 

V/ WWAJ 

P]y^“ q| 

+ i /2L: 


‘ft »m 

'cos” 

•1 P(y - n)l 

P \y - 

-V 


X - i ■ 

1 

= (y - 4 ) 

) 


gu^^(x,. y) = -(y - q) ,f 

X - 

' 1 


-(y ' 

-n) 



Limits of Integrations 

The limits of integration with regard to x for the section 
waveLdrag coefficients and with. regard to y for the total wing 
wave— drag coefficients are discussed. The. u-^Qmponertts caused by 
each of the elementary source lines are zero at all points outside 
of the respective Mach cones. The expressions for the u— integrand 
in equation (1) are therefore evaluated along the section for values' 
of X beginning at the forward boiindary of the Mach cone. This 
integration gives the sectionf-drag— coefficient ccaaponents. In order 
to obtain .the ‘wing wave-drag coefficient, the. s ecti oneway e— drag- 
coefficient components obtained from the rospoctive u-expresslons are 
evaluated along the wing span for values of y contained within the 
Mach cone. Table I refers to one side of the. wing (x and y positive) 
and shows the limits of integration for x and y for the required 
u-expressions for- the wings with sweptbackr and rectangular plan forms 
(fig. 2). 
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... ... . AFEEHDIX B 

F0R4IUIAS. ms WAVE-DMG GdEBnCDfiKTS FOB - UNTi^PFSED : SWEFTBACK 
AND EECTANGULAE WINGS OF BICONVEX PAEABOLIC-ARC 
PHOFILE /il-ZERO LH^ j^pi = cot A = ij 

In. the following analysis the (imntities y and K are. 
employed 'nondimensionally in terms of the semichord , The e'^iiations 
for the drag coefficients in all cases refer to. the real .parts of 
the Indicated expressions* 


Section Wave-Drag Coefficients for Swepfback Wings 

Section vave— drag coefficient without tip effects .— The section 
drag coefficient obtained from eq.uation (i) for the given wing at 
a spahwisa station y .and Mach ntaaber M without the tip effect 
was f oxmd to be as follows : ■ ■ 




3fm! 


J2 


2 cos-1 

Sm’ 


/ 


o(2k; 3 - 3K - 1) cos-1 

. 2m' (K + 1) , 




+ 2)2 _ (Km')2 j> ' 


(Bl) 


where K = ^ and m’ = mp. At the center section, where y or 


m 


K 3 0, eq.uatlon (Bl) beccmeB; 

.2 


32 

00 ^ 3it 


°dco = 




m 


\/m20^ 1 
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IVlien K = 


7 = 


2m 


the forward Mach line frcm the 

> 2m 


— 1 

center section intersects the wing trailing edge, and for y ^ 
equation (Bl) reduces to the following ex^eseion: 


mp — 1 


■ °d«>' - 




m 




Increment in section wave.Hirag coefficient caused hy wing tips .— 
The increment in caused by the tips depends on- the sweep angle/ 
aspect" ratio, and Mach niunher.. The following types occur in an 
unt'apered wing: . * 

1. If A = - 1 -^ — , each tip affects solely its own half of the 

w:lng. In this case the region of the wing affectod is between values 

of y from h ^ — to h. (See fig. 2.) The Increment in section 

mP +1 

wave-drag coefficient at a Mach mmiber M and. spanwise position y 
caused hy the tip was found to he as- follows: 



where the suhscript a Indicates that the x-axis in shifted to the 
tip section, and -v;here y^^’ =; y^P and m‘ = mp. In the plan form 
of“tho wing 


y = y^ + h 


In equation (B2) values for y„ are taken from — -r — ^ — to 0. 

, ■ , ^ mp + 1 

om 

II. If A < — == — , the tip on the opposite wing panel 
mp + 1 

contributes an Increment in in addition to that discussed 



undOT type I. (Sea fig, 2.) The increanent in at a section caused by the opposite 

tip obtained in the foiloulng foim; 




kb 


r2 


12 


(14 + - 3h‘yt’ >. ~ 


cosh' 


_iyb* -2(h» -m»)‘ 


m*y-h' 


3mr3|^i<^ - 1 


+ 3m' 


- 6y-b'®h' syjj'jSh*^ - mi»2) 


r 


cos 


-1 yb*l 

fl- 

+ nP?) - g(h 

* — m* ) 

2m< I 

y 1 _ d’ + m'] 

1 


Cb3)^ 


whore' the subscript b indicates that the x-axls is shifted to the opposite tip . sectfbnj 
and where y^^ = y^p, m’ = mp, and h* = hp. In the plan form, of 'the irt.ng : 

«. , s - 

• ' , * f * ■ • 

7 = y^ - h ' ■ - . ■ ; ' ; ; 

The llml'bs for y-^j to be used in eq.uaticai-(B3) depend on the ■value of the aspect • 
ratio A . Thus 


(a) If A < 


. 2m 


>1, the front Mach line from the opposite tip intersects ■tho 

mp + 1 

■t3*alling edge at a value of y^j, = so that values for y^ in equation (B3) are 

taten fTcm h to • 


< 'Z 


I , ■'fc 
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(b) If A < — SSI — < 1, the front Mach line from the opposite 
+ 1 

tip intersects the tip section, and values fctt* in equation (B3) 

are. taken from, h to 2h. 

In cases (a) and {h) discussed tuader type II, the increment in 
AO(i caused hy the adjacent tip discussed under type I is obtained ' 
at spanwise positions of frcm -h 'co 0. 

The total increment in wave-drag coefficient at a section caused 
by the tips is given by . 


Ac^ = Ac^j + Ac^^^ 

and the total wave-drag coefficient at the section is 


Section Wave-Drag Coefficients for Rectangu,lar Wings 

For the rectangular wing, m’ = mp = «> and the f oll.owlng 
equations result frcm equations .(Dl) bo (E13); 

■ Section wave-drag coefficient without tip effect.— For all 
values of y 


= 


i § 

3 ^ 



(b 4 ) 


This result agrees with Ackeret 's theory for the biconvei parabolic - 
arc profile • 


Increment in ceotj.ori . : f coefficient cauaod by wing tine. - 

^ = p bhe tip effect- results solely frcm the adjacent tip and 


S (tf 
" 3^ \o) 


ya' 


— I2J cosh" 


-1 2 


ye^l 




~ 7a’ ^ 

— 2 C08~^ -. ) 

2 1 

- 

J 


(B 5 ) 
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where 


y=' = = (y - ii)3 


In equation (B5)^ values for y are taten frcm 


a 


to O'. 


If A < .the tip on the opposite wing panel contributes an 

increment in c^j^ in addition to Acg^j . The increment in Ac^j^ at a 
section caused by the opposite tip is 


Ac- 


^II 3n3 



+ 2i 


(B6) 


where 


y-jj’ “ y^jP « (y + ii)P 


The values for y-j^’ or y which are req.uired in equation (b 6) for 

O 

A<- are as follows (see fig. 2(b)); 


(a) 

If 

A>1, 

y 

varies 

frcsE 

0 

to 


(t) 

If 

A<i, 

7 

varies 

from 

0 

to 

h 


_ , Wing Wave— Drag Coefficients for .Sweptback Wings 

The integrations of equations (Bl) to (B3) for the 'section wave- 
drag coefficients between the appropriate limits for y (fig. 2(a) and 
table 1) yield the following results for the sweptback wing. 

^^ing wave— drag coefficient without tip effe ct.— The wing wave- 
drag coefficient without tip effect is ~ - 
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3 A' cosh 


-1 A\ + am» - 
A’ls* 


— '/ 4 m* (A* + m* ) -^ — 1) 


3 


2m' ^ aoa' 




3m ' 3 i/itt'2 .1 1 

..(a«'3 . 3A'm'2-A-3) OOB-J- 

2m' ( A • -I- m’ ) 


(B7) 


where A* = AP and m’ »= mP. 


I ncrement In wing wave-drag c oefficient caused by tips .— If 
^ 2m 

A = — — ^ — j the contribution of the tlpe to the wing wave-drag 

mp + 1 

coeffloient is zero. 

If A < — §2; — ■> i then 

8 /'t^2 „ J 1 


AC. 


D « 


_ 2 cos-i 

2ia’ 


,t2 


(^2m^3 _ A', -r A'^) A' jrL'^ + 1) - 2m ' 


.'3 


A*2 


12nL>3 


2 m,’(n’ — A * ) 

?A‘ cosh~l +\jkm'{m' - A') - A' 2 (m'2 _ 3 ,) 

m'A* ^ 


I 

JJ' 


(B8) 


If A < then 
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Aa 


. 8 f [„ 2^,3 oos-i 

\3^^3f72-ri !_ 2m»2 

+ (2m«3 » 3mt2 ^ ^,3) qqs-I ,±: 1? ..-..ggL 

2syf (m* ~ A’ ) 


+ (A’ — 2m’ )(m’ + A’ )^ cos”^ S:!^ ! j: 

a’ + A' 

— (A’ + 2m’)(m’ — A' )2 cos""^ ^^^ L.. 7 ' ... ^ 

m’ — A' 




12m’- 


6qi' 2 _ ^s2 ^ 


coali"^ i — 4m’ 2 \fx — A ’2 
A’ ’ 


+ 3A’2 cosh. ^ 2^’-* ” , A’ _ ^ A’\/4m’(m’ — A’) — a* 2 (m.’2 _. ^ ) (B9) 

m’A.’ I / 


Total vlng vave-drag coefficient .— The total i/ing wave-drag 
coeffloient is obtained as the sum^ 


where the oanponents Cp and ACj) are calculated from the foregoing 

equations for the wing wave— drag coefficient appropriate to the aspect 
ratio of the wing. -i; ■ ■ ' 


Wing Wave— Drag Coefficients for Rectangular* Wings 

For the rectangular wingj m’ = In^ = <» . The following results 
for the wave— drag coefficient of the rectangular wing aiay bo obtained 
either by integrations of equations (B4) to (B6) between the appropri- 
ate limits for y (see fig. 2(b)) or by substitution of m =a> in 
equations (B7) to (B9); 

Wing wave-drag coefficient without tip offect .— The wing wave— drag 
coefficient without tip effect is 

c. =i6m2 
3P\c/ 


(BIO) 
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This result agrees with Ackeret’s theory for a hiconvex parabolic- 
arc profile . 


Increment in ving wave-drag cooff icient caused by tips The 
increanent in contributed by the tips is obtained in the foiiii 



(BU) 


where 


A« * AB< 1 


If 



6 q.uation. (B3JL) shows that = 0. 
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Limits of Integration 



®Tho ocnpononte and «« aero if A - j ; thoroforo, the integration 


liiuita for these u-oosgionents refer to oaaea In vhich A < ^ ^ Y‘ 

^he oanponents and ^ ^ are aero if A therefore, the integration lial 

* 2 
those u-caoBonenta refer to cases in which A ^ . 
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Indicates Mach hm 


Figure 1.- Sym^bols for sweptback wings. 





(b) Reef angular w mg I Mach number, les, ^ , 

^ -Indicates Mach l/ne 

Fi^e 2.- CoDfigiirations showing typical Mach lines behind wing leading edge for various 
aspect ratios. Biconvez parabolic-arc profile; no taper. 

CO 

“0 


NACA TN No. 1449 


Secf/on drag coefficieni^ 
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Figure 3.- Typical distributions of section wave -drag coefficients 
along wing span for sweptback wing. Biconvex parabolic-arc 

profile at zero lift; no taper; — = 0.10; a = 45°. 
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B/sfance from cenferj Yj semt chords 

(a) Seqiion-drqg-^oeffic/eni.djsir/bui/on wdhou/ 
Tip effect for several Mach numbers. 



Distance from cenfer^ semichords 

ib) Secfion-draq-cgeff icienf distribution with tip effect : 
Mach number^ 1.2.3. 


Figure 4.- T37pical distributions of section wave-drag coefficients 
along wing span for rectangular wing. Biconvex parabolic-arc 

L = 0.10; A = 0°. 
c 


profile at zero lift; 




(a) 


0 ^ cot A 



1 


^ 0.664. 


Figure 5.- Generalized curves for determining variation of wing wave-drag coefficient with 
aspect ratio for constant sweepback angles, Mach numbers, and thickness ratios. Biconvex 
parabolic -arc profile at zero lift; no taper. 
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(c) 1.00 = cot A -Jyr - 1 = 1.50. 
Figure 5.- Continued. 



(K^pecf-rafto pammefer ^ A fan A 


1.5 ^ cot A 


- 1 = 7.0 


Figure 5.- Concluded. 
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Figure 6.- Generalized curves for determining variation of wing wave-drag.. 
coefficient with Mach number for constant aspect ratios, sweepback 
angles, and thickness ratios. Biconvex parabolic -arc profile at zero 
lift; no taper. 
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Figure 6.- Continued 



Drag- coefficient parameter 
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Figure 6.- Concluded. 
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Figiare 7.- Generalized curves for deterniining variation of wing wave-drag 
coefficient with aspect ratio for constant Mach numbers and thickness 
ratios for rectangular wing. Biconvex parabolic-arc profile at zero lift, 

C 

(Arrows indicate that values for remain constant for all 

100 (t/c)^ 


aspect ratios to infinity.) 
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Figure 8.- Concluded. 
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Drag coefficient^ 



Mach number^ M 
(a) 5 fenderness rafio^ 2.5. 


Figure 9,- Variation of wing wave-drag coefficient with Mach number for different sweepback 
angles with constant slenderness ratios. Biconvex parabolic -arc profile at zero lift; no 

taper; A = O. 2 O 7 cos^A; 7 = 0.10 cos A ; constant wing area. 
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Figure 10.- Variation of wing wave-drag coefficient with Mach number for different sweepback 
angles with constant aspect ratios. Biconvex parabolic-arc profile at zero lift; no taper; 

— = 0.10; constant wing area, 
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Mach number^ M 
(b) /\specf rafto^ 5. 
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Figure 10.- Concluded. 
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